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ABSTRACT 

A nonlinear control system sometimes oscillate at a 
frequency that is an integral submultiple of the driving 
frequency in response to a sinusoidal input. Such a re- 
sponse is undesirable. An analog computer study is under- 
taken to investigate the occurrence of subharmonic oscillation 
in second order systems with nonlinearities characterized by 
saturation and hysteresis effects. 
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I. INTRODUCTION 



A. NONLINEAR SYSTEMS 

A system is said to be linear if its performance obeys 
the "Principle of Superposition". The "Principle of Super- 
position" states that if C^(t) and C 2 C t ) are the responses 
of a system when it is separately subjected to two inputs 
R-^(t) and R2(t), respectively, then for all a and b, the 
response of the system to the input aR^ (t) +bR2 (t) is 
aC-^ (t) + bC2 (t) and this must hold for all inputs. In other 
words, once the response to one type of signal is known , 
the response to any other signal can be deduced directly. 

On the other hand, any system is said to be nonlinear 
if the "Principle of Superposition" does not hold. The 
term nonlinear systems applies to varied types of systems 
which have practically nothing in common with one another. 
As a consequence, there is no way of generalizing the re- 
sponses of nonlinear systems for a class of inputs to the 
response for any other inputs . The response to one type 
of input does not necessarily contain sufficient informa- 
tion about the response to some other input. Even two 
inputs which are of the same waveshape but different in 
amplitude can give responses that are entirely different. 
This makes the study of nonlinear systems quite difficult. 
Even with this difficulty, there is little choice except 
to tackle the problem in some way because nonlinear sys- 
tems occupy a very important role in practical systems. 
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There is a great variety of possible nonlinear char- 
acteristics, and there are different ways by which non- 
linear systems are classified. One classification is of a 
mathematical kind based on the equations involved, however 
the study of nonlinear systems on a mathematical basis is 
extremely difficult. In this study, types of nonlinearities 
encountered in servo systems are considered, namely satura- 
tion and hysteresis effects. The consideration of nonlin- 
ear elements in the operation of control systems have 
practical significance because it takes linear theory nearer 
to reality. 

B. SUBHARMONIC OSCILLATIONS 

A linear system when subjected with a sinusoidal in- 
put produces a steady-state output which is also a sinusoid 
and oscillates at a frequency equal to the input frequency 
but may differ from the input in amplitude and phase. This 
is true for all frequencies. However, for a system with 
suitable nonlinear characteristic when subjected to a sin- 
usoidal input, there is no guarantee that the steady-state 
output is also a sinusoid. In the usual case, the output 
wave contains frequency components which are integral mul- 
tiples of the input frequency. Also, the system can pro- 
duce a response which oscillates at a frequency which is 
an integral submultiple of the input frequency. This fre- 
quency response phenomenon is known as "subharmonic reso- 
nance". When the frequency of the output is 1/n times the 
driving frequency, where n is an integer, then the output 
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oscillation may be called the 1/n order subharmonic oscil- 
lation. When the output oscillates with the input fre- 
quency, the output oscillation is called the fundamental 
oscillation. 

The subharmonic. phenomenon is of a resonant nature and 
because of this, the subharmonic component is often of large 
amplitude which causes the system output to bear little re- 
semblance to a sinusoid of the input frequency. The output 
is considerably different from that obtained when the sub- 
harmonic phenomenon does not occur. For systems in which 
the output is intended to reproduce the input as closely as 
possible, as in the motion of a loudspeaker cone or of a 
servo motor, this phenomenon can cause a distortion that 
cannot be tolerated. In the field of feedback control sys- 
tems, it is interesting to know the conditions under which 
subharmonic oscillations may occur. Subharmonic oscilla- 
tions can give false information which can be troublesome, 
so it is desirable to eliminate the possibility of the 
occurrence of subharmonic oscillations. This is true in 
the case of recording seromechanisms which are frequently 
subjected to periodic inputs. The possibility of the occur- 
rence of the subharmonic oscillations reduces the reliability 
of the systems. 
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II. DUAL- INPUT DESCRIBING FUNCTION 



The requirement that a nonlinearity be described by 
an equation in the frequency domain and that this equation 
be compatible for use with the transfer function of the 
linear components gave rise to approximating the effects 
of the nonlinear component by using a linear "approximate 
transfer function" or describing function. The approxima- 
tion is made by defining the describing function in terms 
of the Fourier series for the component response to a sin- 
usoidal input. 

The describing function is defined as the ratio of 
the magnitude of the fundamental term in the Fourier se- 
ries for the output wave to the magnitude of the input 
sinusoid at a phase angle which is the angle between the 
two sine waves and for all permissible amplitude and fre- 
quency of the input wave [Ref. 9]. In mathematical form, 
let the input be ASintot and the fundamental frequency term 
in the Fourier series for the output waveform is 

F (A ,to) Sin [tot + (A, to)], 

then the describing function for the nonlinear component 
is 



DF = 



/4>CA,to) 
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Therefore, the describing function is but a mathemati- 
cal linearization of a nonlinearity in the presence of a 
sinusoidal input. In analysis this sinusoid is physically 
identified as the system limit cycle and the describing 
function is used to predict its amplitude and frequency. 

The use of the describing function is based on the assump- 
tion that harmonics generated by the nonlinear characteris- 
tics are sufficiently filtered before being fed as input to 
the nonlinear element. 

An extension of the describing function technique is 
the consideration of a more complex signal fed to the non- 
linear element. The first extension is the use of dual- 
frequency signal, 

= a cos(wt+<}>) + b cos nut 

as input signal to the nonlinearity with the restriction 
that a, b, and n must be real. The dual-input describing 
function can be defined as the ratio of the amplitude of 
the desired frequency component in the output waveform to 
the amplitude of the component of the same frequency in 
the input waveform. Therefore, in the determination of 
the dual-input describing function it is necessary to find 
the output waveform of the nonlinear element for an input 
composed of two sinusoidal signals and analyze the output 
to find the components of a particular frequency. In the 
systems considered here, these signals are composed of 
that component due to the input signal at a fixed frequency 
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and the component due to the fed back signal which may be 
of different frequency from the input signal. It is again 
assumed that the harmonics of the signal fed back are ade- 
quately filtered. The concept of the describing function 
is thus extended to allow for an input signal composed of 
two sinusoidal signals of different amplitudes and fre- 
quencies . 

To illustrate the derivation of the dual-input des- 
cribing function, consider the cubic characteristic, 

y = x 3 

Let the input be of the form 

x = A sin mot + B sin (iot + 0) 



where 



A and B are independent amplitudes, 

(o and nco are any two values of frequencies , and 
0 is an independent phase angle. 

Substituting the value of x in the nonlinear equation, the 
result is, 

y = [A sin niot+B sin (iot + 0)] 3 

= A 3 sin 3 ncot + 3A 2 Bsin 2 mot sin(cot+0) 

+ 3AB 2 sin mot sin 2 (iot + 0) +B 3 sin 3 (oot + 0) 
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Using trigonometric identities to express the above ex- 
pression as a series of sine waves of fundamental and 
harmonic frequencies results in 



y = A 3 



[I sin 



nu>t - sin 3ntot 



+ 3A 2 B 



j sin(tot+0) - j sin{ (2n+l) tot+0} 



+ j sin{ (2n- 1) tot- 0 } 



+ 3AB 2 



1 1 

2 sin ntot - sin { (n+2) tot+20 } 



sin { (n- 2) tot- 20 } 



+ B 



3 1 

sin (tot + 0) - ^ sin 3(tot + 0) 



or 



y = (A 2 + 2B 2 ) sin ntot - sin 3ntot 



+ ^ (B 2 +2A 2 ) sin (tot + 0) •- sin 3(tot+0) 



3A : 



sin{ (2n- 1) tot- 0 } - sin{ (2n+l) tot + 0 } 



3AB : 



sin{ (n+2) tot+20 } + sin{ (n- 2) tot- 20 } 
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By inspection of the above expression, the sum of the first 
and third terms is the fundamental input function equation, 
x multiplied by the factors ^-(A 2 + 2B 2 ) and j(B 2 +2A 2 ) . For 
the other terms, there are values of n which have to be 
considered as will be explained in succeeding paragraphs. 
Consider the following cases: 

(1). Non-harmonically related input sinusoids. In 
this case n is irrational, i.e., n is not expressible as 
a ratio of two integers. In the problem being considered, 
n / 1/3, 1, or 3. For the frequency na)/ 27 T, the dual-input 
describing function is : 



N a (A,B) 



j A(A 2 + 2B 2 )sin noot 
A sin nwt 



N a (A,B) = | (A 2+ 2B 2 ) 

and for frequency w/2tt, the dual-input describing function 
is : 

j B(B 2 + 2A 2 ) sin(u)t+0) 

N g (A,B) = B sin(a)t + 0) 

Ng(A,B) = |(B 2 + 2A 2 ) 

The resulting dual-input describing functions, N A (A,B) and 
Ng(A,B) are both independent of the phase angle, 9 but 
dependent on the input signal amplitudes A and B. 

The dual-input describing function is modified if n 
assumes a value such that the other terms in the series 



17 



have the frequency of either of the input signals co/ 2 tt or 
nw/27r, i.e., n = 1/3, 1, or 3. 

(2). Harmonically related input sinusoids. 

(a) n = 3. For this case there are two ad- 
ditional terms in the output waveform, y that is considered, 
i.e., the fourth and eight terms, such that the dual-input 
describing function becomes 



N 



a (A,B,3,6) = 



(A 2 +2B 2 ) sin 3wt 



R 2 

— sin (3cot + 30) 



A sin 3wt 



The above expression shows that varied phase angles 
are associated with the terms, so the phasor representation 
is utilized, then 

f A (A 2 + 2B 2 ) e 0 - |1 e jt 3 “ t + 30 ’ 
N a (A,B, 3,0) = i 4 - 



N A (A,B,3,e) = | (A 2 +2B 2 ) - 21 e 030 
Similarly for 

22 (B 2 *2A 2 ie 3 f“ t+e ’ - |22i e jt wt - 20:) 

N (A,B,3,e) = 2 2 

Be j (« t+e ) 

N B (A,B,3,e) - | (B 2 +2A 2 ) - 222 e' j30 

The resulting dual-input describing functions are complex 
and dependent on A, B, and 0. The resulting DIDFs for the 
other two cases, n=l and n=l/3 as shown in the succeeding 
derivations are also complex and dependent on A, B, and 0. 
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(b) n=l. For this case, and combine to 
give terms of frequency o)/2tt. Proceeding in the same 
manner as for the case n=3, 



N a (A,B,1,0) = 



J A(A 2 +2B 2 )e jwt + p. (B 2 + 2A 2 )e j(a)t + 0:) 



Ae 



jcot 



3A 2 B j (cot - 0) , 3AB 2 j (a)t+20) 



Ae 



j wt 



3 ^* 2 x^ 2 ^ + 4®. 2 rB 2 + 2A 2 'le J0 + e‘j 0 



N a (A,B, 1,0) = ^ (A Z + 2B 2 ) + ^ (B 2 + 2A 2 ) e J + 

+ Mi e ^' 20 

4 e 



N b (A,B,1,6) 



| A(A 2 + 2B 2 )e :it0t + (B 2 + 2a2) e j (^t + 0) 

Be ^ wt + 0 ) 



3A 2 B j (w t ~ 0 ) 3AB 2 j(wt + 20) 

4 e 4 e 

BeJ (wt + 03 



N b (A,B, 1,0) = M (A 2 + 2B 2 ) ' e"** 0 + |(B 2 + 2A 2 ) + e~ j26 

+ m e j 0 



(c) n = 1/3. This case is identical with the 
case for n=3. 
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N a (A,B,1/3,0) 



^ (A 2 + 2B 2 )e j(wt/3: > 1 
Ae j (wt/3) 



N A (A,B, 1/3,0) = | (A 2 +2B 2 ) - ^ 9 

-Al e j«t + 3B (B 2 +2a2) 
N b (A,B, 1/3,0) = Be j (ut+0) 

N b (A,B, 1/3, ) = - ^ e' j0 + | (B 2 +2A 2 



3A 2 B (wt/ 3+0) 

~T 6 



e j (wt+0) 



) 
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III. SYSTEM NONLINEARITIES 



In most practical feedback systems, the nonlinear 
element is in the forward path. The basic system con- 
sidered in this study is shown in Figure Ill-la, where N 
is the nonlinear element whose gain depends on the magni- 
tude of the input to this element and G(s) is the linear 
frequency dependent portion. 

A. SATURATION 

The nonlinearities considered here are saturation, 
relay with hysteresis and backlash. Most physical systems 
contain certain nonlinearities. The most common nonlinear 
phenomenon is characterized by a saturating or limiting 
characteristic. Limiting is almost universally present in 
control systems since most instrumented signals can take 
values only in a bounded range. Many error detectors, 
such as a resolver or synchro differential have a restrict 
ed range of linearity. The saturation phenomenon is also 
present in amplifiers, whether electron tubes, transistors 
magnetic amplifiers or other devices are used. The input- 
output characteristic is shown in Figure Ill-lb. 

B. NON- IDEAL RELAY 

The other nonlinearity is the relay with hysteresis 
which is characterized by the curve shown in Figure III-lc 
The hysteresis delays the switching operation until the 
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NONLINEAR ELEMENT LINEAR ELEMENT 






(c) (d) 



NONLINEAR LINEAR NONLINEAR 




(e) Basic system with two nonlinear elements. 



FIGURE 1 1 1 - 1 : System block diagrams. 
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output is past the desired zero point. A system using this 
type of relay must have a limit cycle in steady state. 

The relay is a device to apply the power which drives 
the load. In a relay control system, the power amplifier 
is a relay device. The relay amplifier is desirable be- 
cause it can be simple, rugged, compact, and relatively 
cheap, while meeting high load power requirements. Ex- 
amples where relay control systems are typically applied 
are in space vehicle attitude control and aircraft and 
missile adaptive control systems. 

C. BACKLASH 

Another commonly encountered nonlinear phenomenon in 
most mechanical linkages is backlash. In control systems 
these are present in most mechanical gears. This is true 
even for new ones and with use, the backlash increases. 

In simple physical terms, backlash in gears is the free 
space between adjacent teeth. It occurs whenever there 
is inexact linkage between mechanical parts, i.e., when 
there is slack between mating gears. The backlash char- 
acteristic is shown in Figure Ill-Id, where A represents 
the amount of backlash. 

In one way or another, linear and nonlinear systems 
will encounter the backlash problem if they have a gear 
train. In the field of control systems, a system is in- 
cluded in the class of unstable systems if it continuously 
oscillates with constant amplitude and zero input (limit 
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cycle). Thaler and Pastel [Ref. 9] show that for a second 
order servo with backlash, the existence of the limit cycle 
depends on the value of the damping factor, For z, 0.29, 

the system is stable so that for z, < 0.29 the effect of 
backlash is to cause an unstable character. If the phase 
portrait is drawn for the system, the phase portrait will 
show the backlash effect near the origin. For systems 
where it is desired to have the phase trajectories termin- 
ate at the origin, the backlash effect causes serious 
trouble . 

The three nonlinearities discussed above were simula- 
ted on the electronic analog computer. The operation of 
the simulators for the nonlinearities are explained in 
Appendix B. The characteristics of the nonlinearities were 
produced very nicely as shown in Figure III- 2. 

Figures 1 1 1 - 3 , III-4, and 1 1 1 - 5 show the block dia- 
grams of the second order systems with the indicated non- 
linearity which were simulated on the analog computer. 
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FIGURE 




(a) Saturation 





1 1 1 - 2 : Simulation output for nonlinearities 
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FIGURE III-3a: Block diagram of second order system with 

saturation . 




FIGURE III-3b: Analog simulation using digital logic for 

simulation of saturation. 



26 




FIGURE III-4a: Second order system with non-ideal relay. 




FIGURE 1 1 1 - 4b : Analog simulation of system with relay. 
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FIGURE III-5a: Second order system with, two nonlinearities. 
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FIGURE III-5b: Analog simulation of system with saturation and backlash. 



IV. SUBHARMONIC GENERATION 



Consider the basic system shown in Figure Ill-la with 
a nonlinear element in the forward path. If the nonlinear 
element is replaced by an element of unity gain, the result 
ing linear system is assumed to be stable. 

As stated earlier, the "Principle of Superposition" 
applies to linear systems. The response of the closed loop 
system to a given signal is independent of any other signal 
which may be applied at the same time. Mien a nonlinear 
element is included in the feedback loop this is no longer 
the case. The response of the nonlinear element to a given 
signal is modified if another signal is applied at the same 
time. The effect of the other signal may be considered as 
varying the gain of the nonlinear element for the primary 
input signal. The case of interest is when the additional 
signal input is also sinusoidal with a frequency that is 
an integral multiple of the primary input frequency. In 
In this case, the output component of the angular frequency 
a) may be modified in amplitude and phase by the presence of 
the additional input. This modifies the loop gain of the 
system. The phase change is caused by the intermodulation 
products generated by the nonlinear element. 

The stability of the closed loop system is affected by 
varying the gain of the nonlinear element since the nonlin- 
ear element is a part of the feedback loop. For a range of 
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sinusoidal inputs, the response of the nonlinear element 
may be so modified that the closed loop is rendered oscil- 
latory. Where the system is such that the open-loop fre- 
quency response plotted on the complex plane do not 
intersect with the negative real axis, the nonlinear ele- 
ments have to introduce an adequate phase change in order 
that stable oscillation can occur. This implies that the 
frequency of the oscillation must be related to the input 
frequency. What is of interest here is where the fre- 
quency of the oscillation is an integral submultiple of 
the primary input frequency in which case subharmonic os- 
cillation is generated. 
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V. DISCUSSION 



The systems shown in Figures 1 1 1 - 3 , III-4, and 1 1 1 - 5 
with the indicated nonlinearities were simulated to in- 
vestigate the occurrence of subharmonic osciallations and 
to observe the effect of changing the variable parameters 
in the systems such as, the damping factor, £ and forward 
linear gain, K m in Figures 1 1 1 - 3 and 1 1 1 - 5 and the amount 
of hysteresis, 2c in Figure III-4, in the generation of 
subharmonic oscillation. 

A. SYSTEM WITH SATURATION NONLINEARITY 

Figure III- 3a shows a second order system with a sat- 
uration nonlinearity. The nonlinear characteristic is 
defined as : 

e = e . 2 < e - <2 

0 1 1 

e = 2 e . > 2 

o 1 

e=-2 e - < - 2 

o 1 

and the gain of the linear portion of the saturation is 
unity. The transfer function of the linear part is, 

G(s) = s(s T-074T 

Considering the linear part, the closed loop transfer 
function is 

K m 

s" 2 + 0.4s + K 
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with characteristic equation 

s 2 +0.4s+K =0 

m 

The natural frequency, and the damping factor, x> is 
determined from the following relationships: 



w 



n 



= a r 

m 



C 



0.4 

2 “n 




2 



In the simulation of the system in the electronic analog 
computer, was varied from 1 to 100 giving a variation 
in c from 0.2 to 0.02. Typical results obtained by the 
simulation are shown in Figures C-l through C-5. 

Ogata [Ref. 6] discussed an analysis of subharmonic 
responses in systems with a nonlinear element whose input- 
output characteristic curves shows saturating or limiting 
characteristics and showed that for nonlinear characteris- 
tic curves that are skew symmetric, odd-order subharmonic 
oscillations can occur. 

The "28-degree criterion" for predicting responses of 
saturating feedback systems was introduced by Douce [Ref. 1] . 
In the complex plane, two lines are drawn from the point 
s = -1 on each side of the negative real axis and each line 
makes an angle of 28°. This is illustrated in Figure V-l. 

The region bounded by the two lines represent the critical 
region for a system with saturation nonlinearity. Utilizing 
this, the Nyquist diagram of the linear part of the system 
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Ttvi 




FIGURE V-l: The 28-Criterion (Douce [Ref. 1], Shaded area 
shows region where 1/3 subharmonic can be gen- 
erated in systems containing a limiter. 

Dotted lines represent open loop frequency 
response of linear system with indicated value 

of K . 
m 
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provides sufficient information to predict possible occur- 
rence of subharmonic oscillations. Subharmonic resonance 
is possible when the open loop frequency response of the 
linear portion enters the critical region. Figure V-l 
shows the open loop frequency response of the linear part 
for different values of K . Also indicated is the damping 
factor. The figure shows that for K m < 4.S (£> 0.0944), 
the open loop frequency response of the linear portion 
lies outside the critical region so that subharmonic os- 
cillation cannot occur. For K > 4.5, for example K = 50, 
the linear open loop frequency response is inside the crit- 
ical region. The results of the analog computer study 
agree accurately with the "28-degree criterion" as can be 
seen in Figure V-2 which shows the domain of 1/3 subhar- 
monic oscillation as a function of input frequency and 
damping factor, £. The figure also shows the effect of 
input amplitude, i.e., for greater input amplitude, the 
domain decreases and the system has to be more lightly 
damped to sustain subharmonic oscillation. 

Subharmonic oscillation of 1/3 order was the only one 
observed for the system with skew symmetric saturation non- 
linearity. Figure V-3 shows the domain of 1/3 subharmonic 
as a function of input amplitude and input frequency. It 
can be seen that the range of input frequency where sub- 
harmonic oscillation is sustained decreases with increasing 
input amplitude. The same observation can be seen in 
Figure V-2. It can be deduced that subharmonic oscillation 
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INPUT FREQUENCY, rad/sec 

FIGURE V-2: Area for occurrence of 1/3 subharmonic for symmetric saturation. Grap 

shows relationship between input frequency and damping factor and inpu 
amplitude. Darkened marks are obtained by experiment. 



INPUT AMPLITUDE 3 volts 6 volts 9 volts 16 volts 
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Table V-la: Range of input frequencies that sustain subharmonic oscillation. Data 

obtained from results of analog computer simulation. 



n Input Frequency Input Frequency Input Frequenc 

rad/sec Range, rad/sec Range, rad/sec Range, rad/sec 
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Table V-2 : Output amplitude of system with saturation 

X, = 0.02. 
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can occur only for a range of input amplitude. For a small 
input amplitude such that the error signal to the nonlinear 
element is less than the saturation level, no subharmonic 
oscillation can possibly occur because the system will just 
be linear. Also shown in the figure is the variation of 
output amplitude. For a particular damping factor, the 
example l , = 0.02 in Figure V-4, the output amplitude varies 
as shown by the curve irregardless of the input amplitude. 
The output amplitude is dependent on the range of input 
frequency that the subharmonic oscillation is "locked on". 
Figures V-4 and V-5 show the output amplitude variation 
for different values of C . Ogata [Ref. 6] concluded that 
subharmonic oscillations do not occur as response of the 
output of the nonlinear feedback control system when the 
output- input amplitude ratio for a given amplitude of a 
sinusoidal input is 2 'v 3 db or less. This is equivalent 
to a ratio of 1.256 1.41. The results of the simulation 

study do not agree with this criterion as can be seen in 
Figures V-4 and V-5 where the lowest output- input amplitude 
ratios are 1.25 and 1.125, respectively. The possibility 
of subharmonic oscillation is greater for greater output- 
input ratio. 

Subharmonic oscillation of 1/2 order was observed when 
the saturation nonlinearity is skew asymmetric. The sat- 
uration characteristic is 

e = e . -2<e.<4 

o 1 1 

e o = ‘ 2 e i < - 2 

e = 4 e . > 4 

o 1 
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Output amplitude, volts 



60 




1 Z "3 4 5 6 7 6 9 To 

Input Frequency, rad/sec 

FIGURE V-5: Output amplitude for 1/3 subharmonic due to 

saturation. Input amplitude = 16 volts. 
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The domain of 1/2 subharmonic is shown in Figure V-6. Fig- 
ures C-l through C-5 show the signal waveforms for 1/2 and 
1/3 orders subharmonic oscillations and also the waveforms 
of the error signals. Figure C-2 shows the build up of 
1/2 order subharmonic and Figure C-4 shows the decay of 
1/3 order subharmonic to fundamental oscillation. 

West, Douce and Livesly [Ref. 12] presented a method 
of using the dual-input, describing function technique to 
study subharmonic oscillation. Figure V-7 is the describ- 
ing function loci which shows the response of a saturation 
nonlinearity to a sinusoidal signal of amplitude B and of 
a signal at three times this frequency and of amplitude A. 
Ref. 12 tabulates other describing function loci for dif- 
ferent values of amplitude B. 

A subharmonic oscillation of 1/n order where n is an 
integer will oscillate at a frequency that is 1/nth of the 
input frequency. To investigate this subharmonic, let the 
input to the nonlinearity be 

A cos nwt + B cos (wt+0) 

The signal A cos nwt is assumed to be known and the des- 
cribing function for the signal B cos (wt+0) is superimposed 
on the open loop frequency response locus of the linear 
part. The value of B, such that the describing function 
cuts the frequency response locus for a frequency, w is 
the amplitude that corresponds to an overall open loop gain 
of unity. So this gives a possible steady-state amplitude 
of the resultant subharmonic oscillation. 
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FIGURE V-6: Domains of subharmonic for a symmetrical saturation input ampli- 

tude = 6 volts. 



The amplitude and frequency of possible 1/3 order 
subharmonic oscillation is found as follows. Assume that 
the amplitude, A of the input signal is known. For the 
system to oscillate continuously at a frequency equal to 
one- third of the input frequency, the overall open loop 
gain must be unity. Then the conditions for continuous 
oscillations may be evaluated by means of the describing 
function technique. The input to the nonlinear element 
is 

A cos 3wt + B cos (wt+0) 

Know the value of A, the describing function loci that 
correspond to this value of A are superimposed on the open 
loop frequency response locus of the linear system. Refer- 
ing to Figures V-7a and V-7b this shows the describing 
function loci for A = 4 and A = 8, respectively. The open 
loop frequency response locus of the linear portion with 
a transfer function of 

Gl » = j ^ j - Mvtrrr 

is superimposed. The point at which they pass through the 
point representing the frequency response at the subhar- 
monic frequency determines the value of B, which is the 
amplitude of the subharmonic oscillation. 

Figure V- 8 shows the output amplitude for 1/3 sub- 
harmonic oscillation for a second order system with sat- 
uration nonlinearity where the saturation level is ± 1. 

The curves were obtained from the results of the analog 
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Q-(iu-A ^ 

jco(jto+0.4) 




FIGURE V-7a: Describing function loci for A = 4. 

Shows response of saturation nonlin- 
earity to a sinusoidal signal of 
amplitude B and of signal AE three 
times this frequency and of ampli- 
tude A. [Ref. 12] . 
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FIGURE V- 7b : Describing function loci for 

A = 8. Shows response of sat- 
uration nonlinearity to a sin- 
usoidal signal of amplitude B 
and of signal at three times 
this frequency and of amplitude 
A. Open-loop frequency response 
locus of linear system super- 
imposed. [Ref. 12.] 
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Amplitude of system response and error signal 
for system with saturation level ±1. 
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computer study. Also indicated are points calculated using 
the describing function loci of Figure V-7 and those tabu- 
lated in Ref. 12. 

The difficulty with the dual-input describing function 
technique is the derivation of the dual- input describing 
function itself. If the nonlinear characteristic can be 
expressed as a simple power series or simple output-input 
relationship such as the cubic nonlinearity considered in 
Chapter II, the dual-input describing function can be de- 
termined by simple algebraic analysis. However, in the 
case of practical feedback systems where the nonlinearity 
is rarely expressible as a simple power series, like in the 
case of saturation, relay, and backlash, the algebraic an- 
alysis of the nonlinear element becomes difficult. Graph- 
ical methods are possible but very laborious. Due to the 
complexity of these nonlinearities, the use of computers 
for determining the relevant dual-input describing func- 
tion is necessary. 

For the purpose of this study, the derivation of the 
dual-input describing function for the nonlinearities are 
not included. The DIDF derivation can be treated as a 
separate problem and the results can then be applied to 
any form of linear systems to study the effect of the 
nonlinearity . 

B. SYSTEM WITH NON- IDEAL RELAY 

The other second order system simulated to observe 
the occurrence of subharmonic oscillation contained a 
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double valued nonlinearity. Figure III-4a shows a second 
order system with a relay with hysteresis. Referring to 
this figure, the nonlinearity can be defined by the fol- 
lowing equations. Assume the input to the nonlinearity 
is 



then 



e^ = A sinwt 



e Q = -a 0 < wt < 8 



e = a 
o 



$ < cot < u+3 



where 



o • -1 

$ = sm c 
A 

The conventional describing function [Ref. 9] is 



r = — 
b D irA 



sin ^c 



and the plot of the inverse of this describing function 
for a = 1 (constant) and three different values of c, where 
2c is the amount of hysteresis, is shown in Figure V-9. 
Superimposed on the -1/G^ curve is the plot for G(jco) of 
the linear portion of the system. The figure shows that 
for a particular amount of hysteresis, the point of inter- 
section between the G(ju>) curve and the relevant -l/G^ 
curve represents a limit cycle whose amplitude and fre- 
quency can be determined easily. For example, consider 
c = 0.2, the limit cycle frequency is approximately 1.8 
rad/sec with an amplitude of 0.46. The Bode plot of the 
simple linear portion is shown in Figure V-10. 
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Amplitude, db 




Phase, degrees 



FIGURE V-9 : Describing function for relay with hysteresis. 



52 



-so 




53 



FIGURE V- 10 : Bode plot of G(s) 



The system is Figure III-4a has a sinusoidal input of 
known amplitude and frequency. The output of the relay, 
f(t) is assumed to lose its higher frequency components due 
to the low pass characteristics of G(s) . Then, the input 
to the nonlinearity can be written as 

A cos nwt + B cos (wt + 0) 

where A cos nwt is the primary input signal of known am- 
plitude and frequency. The dual- input describing function 
of subharmonic oscillation developed by West, Douce and 
Livesly [Ref. 12] can be used if the dual-input describing 
function is available. 

In the analog computer study for the system, subhar- 
monic oscillation as high as 1/2 order and as low as 1/13 
order were observed. The subharmonic oscillation once ini- 
tiated is stable and stays "locked on" a reasonable range 
of frequency. There is an overlapping of domains of dif- 
ferent subharmonic order, as shown in Figure V-ll. The 
subharmonic oscillation that is stable is the subharmonic 
order initiated first. It can.be seen from Figure V-ll 
that odd-order subharmonics have a much larger domain than 
even-order subharmonics and that the domain decreases for 
lower odd-order subharmonics. The figure also shows the 
amplitude of the output which is the amplitude of the sub- 
harmonic oscillation. This is the amplitude curve for all 
subharmonic modes. If the output oscillates at a frequency 
of 2.45 rad/sec, then the output amplitude is 0.225 irre- 
gardless of the subharmonic mode. 
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Table V-4: Simulation Data for Non-Ideal Relay; 

a = 1 , c = 0 . 2 
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Table V-5: Simulation Data for Non-Ideal Relay (continuation) 
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Input Amplitude, volts 
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s * + s + 1 



, Boundary points obtained 

l by analog computer. 



» 




1 2 5 ? 5 ^ 

OUTPUT FREQUENCY, rad/sec 



-11 : Domains for subharmonic output for system 

with relay: a = 1, c = 0.2. 
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Figure V-12 shows the domain of subharmonic oscilla- 
tions for a = 1 and c = 0.1 adapted from the work of Sawaragi 
and Akashi [Ref. 8]. The points plotted represent data that 
were- obtained from the analog computer experiment. 

Figures C-6 through C-16 show the simulation results 
for the system with non- ideal relay. In the system with 
saturation nonlinearity, subharmonic generation was not 
self starting. The occurrence of subharmonic oscillation 
was dependent on the initial conditions, such as, the sud- 
den change of input frequency or input amplitude. For the 
system with non- ideal relay nonlinearity, subharmonic oscil- 
lation was self starting. This can be seen in Figure C-6 
which shows the system response changing from fundamental 
oscillation to 1/3 subharmonic oscillation. Figure C-6 
also shows the relationship between the amplitude of the 
output when oscillating at the fundamental frequency and 
when oscillating at the subharmonic frequency. Figure C-10 
shows the change from subharmonic to fundamental oscilla- 
tion, i.e., from 1/5 order subharmonic to fundamental. The 
amplitude of the output at the fundamental frequency is 
practically negligible compared with the output amplitude 
when oscillating at a subharmonic frequency. Figures C-7 
through C-16 show the changing modes of the subharmonic os- 
cillation due to change in input frequency. There were in- 
stances when the output oscillates at two different subharmonic 
modes at the same input frequency. Figure C-16 shows the 
output oscillating at two different subharmonic modes ; name- 
ly, 1/11 order and 1/13 order at the same input frequency. 
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C. SYSTEM WITH TWO NONLINEARITIES 

The third system simulated to investigate the occur- 
rence of subharmonic oscillation contained saturation and 
friction controlled backlash nonlinearities as shown in 
Figure III-5a. It is a two loop system with a tachometer 
feedback for the inner loop to provide the necessary damp- 
ing to make the system stable. 

Using standard procedure, the characteristic equation 
is 



It is clear 
and , i.e 



s 2 + K K.s + K = 
m t m 

that the damping factor, 



C 




0 

C is dependent on K m 



Han [Ref. 4] presented a procedure to predict the 
existence of limit cycles in a system with two nonlinearities 
such as the system shown in Figure 1 1 1- 5a. For the purpose 
of observing the occurrence of subharmonic oscillation in 
the system, assume that the system oscillates with a limit 
cycle so that the input to the first nonlinearity is of the 
form 



A cos mot + B cos (wt+0) 

where A cos mot is the primary input signal of known ampli- 
tude and frequency. The linear frequency dependent ele- 
ment, G(s) is a low-pass filter which is assumed to adequately 
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attenuate the higher frequencies so that harmonics can be 
neglected in It also follows that the input to the 

second nonlinearity is a sinusoid. 

The analog computer simulation was performed for dif- 
ferent values of and K t and the results are tabulated 
in Table V-6 and V-7. Figures 017 through 0 23 show the 
signal waveforms obtained. Figure 017 and 018 show the 
waveforms for 1/3 and 1/5 orders, respectively. It can 
be seen that the signal e^ n which is the input to the back- 
lash nonlinearity oscillates at the subharmonic frequency, 
i.e., in Figure 017, e^ n is oscillating at a frequency 
that is 1/3 of the input frequency and in Figure 018, at 
1/5 of the input frequency. The system output, 0^ oscil- 
lates at the same frequency as e^ n with a phase difference 
of 180°. This shows that the subharmonic oscillation is 
due to the effect of the saturation nonlinearity. The ef- 
fect of the backlash nonlinearity can be analyzed by using 
the conventional describing function since the input to the 
backlash nonlinearity, e^ n is only composed of one sinusoid- 
al signal. Figure 019 shows the change of fundamental os- 
cillation to 1/3 subharmonic oscillation and Figures 020 
and 021 show the build up of subharmonic oscillations. 

It can be seen also from these figures that the subhar- 
monic oscillations are self starting as compared with the 
system with saturation nonlinearity only. 

The domains of subharmonic oscillation as a result of 
the simulation study are shown in Figures V-13, V-14, and 
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V-15. The domain of the 1/3 subharmonic oscillation is of 
the order of the natural frequency. For lower order sub- 



harmonic modes, i.e., 1/4, 1/5, 1/6, 1/7, 
oscillation occurs at a higher input frequency 



subharmonic 
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Table V-6: Domain of subharmonic for system with two 

nonlinearities, £ = 0.283, & n = 14.14 
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3.0 9.0 9.6 21.0 25.5 30.0 

Input Frequency Frequency Frequency Frequency Frequency Frequenc; 

Amplitude Range Range Range Range Range Range 

Order rad/sec rad/sec rad/sec rad/sec rad/sec rad/sec 
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Table V-7: Range of input frequency for subharmonic oscillation (two nonlinearities) 
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FIGURE V-14: Domains of subharmonic for system with saturation and backlash 

t=0 . 2 . w ~ 10.0 
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FIGURE V-15: Domains of Subharmonic oscillations for system with two nonlinearities 

£ = 0.141; co = 14.14 rad/sec 



VI. CONCLUSIONS 



The investigation of subharmonic oscillation in non- 
linear systems was undertaken by means of an analog com- 
puter study on three nonlinear second order feedback 
control systems. The following conclusions are summarized 
from the result of this study. 

When the output of the system with saturation non- 
linearity is oscillating with the input frequency, sub- 
harmonic oscillation cannot occur unless some kind of sudden 
disturbance is given to the input. This sudden disturbance 
may be a sudden disturbance may be a sudden increase in 
input amplitude or frequency. This is not the case for 
the systems with hysteresis nonlinearity. In the system 
with non-ideal relay and backlash characteristics, sub- 
harmonic generation was self starting. 

The possibility of the occurrence of subharmonic os- 
cillation is great when the output- input amplitude ratio 
is large. In the case of saturation nonlinearity, there 
is a minimum output- input amplitude ratio for subharmonic 
oscillation to occur, i.e., the output-input amplitude 
ratio is greater or equal to 1.125. For the system with 
non-ideal relay, subharmonic oscillation can occur for 
output-input amplitude ratio less than unity. 

The "28-degree criterion" introduced by Douce [Ref. 1] 
gives accurate results in predicting the possibility of 
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subharmonic responses in a system with saturation nonlin- 
earity. This is a useful technique in the elimination of 
subharmonic oscillation for systems with a limiting non- 
linearity . 

Subharmonic oscillation once initiated is sustained 
over a reasonable range of frequency. There is an over- 
lapping of these frequency ranges for different subhar- 
monic modes. In the region where two or more orders of 
subharmonic oscillation can occur, the order that occurs 
first dominates and tends to continue and is quite stable. 
The frequency of the subharmonic oscillation has a minimum 
depending on its order. The frequency for 1/3 order sub- 
harmonic is usually in the neighborhood of the natural 
frequency of the system. Subharmonic oscillation is de- 
pendent on the input signal amplitude. It can occur only 
within a certain range of signal amplitude. It can be con- 
cluded that subharmonic oscillation can occur only with a 
certain range of input amplitude and frequency. 

The system with a limiting nonlinearity is dependent 
on the damping factor for subharmonic oscillation to occur. 
The system has to be very lightly damped, i.e., C < 0.1 in 
order for the system response to oscillate at a subharmonic 
frequency. If the input amplitude is increased, the system 
has to be much more lightly damped for subharmonic oscilla- 
tion to occur. 

The dual- input describing function technique developed 
by West, Douce, and Livesly [Ref. 12] appears to be a very 
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useful technique for predicting quantitatively the ampli- 
tude and frequency of subharmonic oscillation. Recommend- 
ed for future study is the derivation of relevant dual-input 
describing functions (two sinusoidal input) and application 
of the results to any linear system with the relevant non- 
linearity. Quantitative analysis of the occurrence of sub- 
harmonic oscillations should be undertaken. 
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APPENDIX A: PROGRAM SYMBOLS 



(1) Servo Set Potentiometer 



( 2 ) 




o 



o 



e = Ke . where K < 1 
o 1 



e 

o 



On * 



10c i2 ) 



(3) Inverter 




(4) Integrator 




e = - 
o 



e . 

i 



e 

o 



( e il+e i2 +1 0 e i 3 ) dt 



(5) Comparator 




Output Logic: 

Algebraic sum of X + Y > 0 

C = -6 volts (true) 

C = 0 volt (false) 

Algebraic sum of X + Y < 0 
C = 0 volt (false) 

C = -6 volts (true) 
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(6) Digital/Analog (D/A) Switch 




Logic Control: 

The switch is closed if either 
-6 volts is applied or the 
input is left open; the switch 
is open if 0 volt is applied. 



(7) OR Gate 



" L 1 + L 2 



(8) AND Gate 



* L 1 + L 2 
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APPENDIX B: ANALOG SIMULATION 



1. Nonlinearities 
(a) Saturation 

Figure B-la shows the saturation characteristic 
curve and Figure B-lb is the circuit diagram for the simu- 
lation using the analog computer. Referring to Fig. B-lb, 
the simulation works as follows : 

(1) e i < - a 

Algebraic sum of Cl < 0 Cl false, CT 
true and also the algebraic sum of C2 < 0 -*■ C2 false, 

C 1 true such that the output of the AND gate = 0 and 
D/Al is open. For Cl true, D/A2 is closed and for C2 
false, D/A3 is open, therefore, 

e = + V, 
o 1 

(2) Algebraic sum of Cl > 0 -*■ Cl true, Cl 
false and the algebraic sum of C2 < 0 C2 false and C2 
true. The output of the AND Gate = -6 volts and D/Al 
is closed. For Cl false, D/A2 is open and C2 false, 

D/A3 is open, therefore, 

e = - me . 
o 1 

(3) 0 < e i < b 

Algebraic sum of Cl > 0 Cl true, Cl 
false and the algebraic sum of C2 < 0 C2 false and C2 
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e, 




FIGURE B-la: Saturation characteristic curve. 




FIGURE B-lb: Analog simulation for saturation. The 

symbol "g" refers to frid of amplifier 
and 100 is the reference voltage of the 
computer being used. 
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true. The output of the AND Gate = -6 volts and D/Al is 
closed. For Cl false, D/A2 is open and C2 false, D/A3 
is open, therefore, 

e o = : me i 
(4) e i > b 

Algebraic sum of Cl > 0 ■+ Cl true, CT 
false and the algebraic sum of C2 > 0 -*■ C2 true and C2" 
false. The output of the AND Gate = 0 and D/Al is open. 

For Cl false, D/A2 is open and for C2 true, D/A3 is 
closed, therefore, 

e » - V, 
o 2 

In the description of operation above, -6 volts represents 
a logical 1 and 0 volts represents a logical 0. The logic 
control for the D/A switches is the output of the compara- 
tor or of the AND Gate . 

(b) Relay with Hysteresis 

Figure B-2a is the characteristic curve of a 
non-ideal relay with hysteresis and Figure B-2b is the 
diagram for analog simulation. The Track-Track (T/T) unit 
is composed of two D/A switches with a common G terminal. 
Figure B-2c is the analog circuit for the T/T unit. Switch 
T is closed and T is open when a logical 1 is applied to 
the input or the input is left open. When the input is a 
logical 0, the switch states are reversed. When the 
switches are open, the resistor switch ends are grounded. 
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FIGURE B-2a: Characteristic curve of relay with hysteresis. 




FIGURE B-2b: Analog simulation of relay with hysteresis. 



R 

T VVW 

R 

T* VWV- 




Q 




FIGURE B-2c: Analog circuit of track-track (T/T) unit. 
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(c) Backlash 

The backlash characteristic curve is shown in 
Figure B-3a and the analog simulation circuit diagram is 
shown in Figure B-3b. The backlash characteristic is 
achieved as follows: 

As e^ increases from zero nothing happens to e Q 
until (e^ - e ) = e' > a. At that point the integrator 
output will begin to change at the rate m(e^ - e ) , and 
e Q will approach (e^ - a) through the feedback correction 
of the loop. As soon as e^ changes direction, (e^ - e ) 
becomes less than a resulting in e" going to zero and the 
output of the integrator holding its value. As e^ con- 
tinues to decrease, e Q remains unchanged until e' becomes 
less than b, at which time, the output of the integrator 
begins to change at the rate m(e^ - e Q ) . Again when e^ 
changes direction, e 1 ' goes to zero and the output of the 
integrator remains at some steady voltage. 

The results of the simulation for the three nonlin- 
earities considered are shown in Figure 1 1 1 - 2 . 

2. Scaling 

When simulating problems on an analog computer, scaling 
is quite important. There are many approaches to this 
technique. What is of interest here is to point out that 
it is necessary to perform magnitude and frequency scal- 
ing on the analog simulation of the nonlinearities simula- 
ted in this study. This is true for the two nonlinearities 
with hysteresis. The simulation for the saturation nonlin- 
earity only requires amplitude scaling. 
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FIGURE B-3a: Backlash characteristic curve. 




FIGURE B-3b: Analog simulation for backlash. 
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FIGURE C-l: 1/2 subharmonic oscillation due to asymmetric 

saturation input = 6 Sin 1.6Trt; = 20 
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FIGURE C-2: Fundamental to 1/2 subharmonic due to asymmetric 

saturation, smm/sec. 
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FIGURE C - 3 : 



1/3 subharmonic oscillation due to symmetrical 

saturation input = 3 Sin 5.6Trt, K = 100 
r m 
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FIGURE C-4: 1/3 subharmonic oscillation to fundamental os- 

cillation input frequency from 2.9hz to 3.0hz 

K =70. 
m 
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FIGURE C-5: 1/3 subharmonic oscillation due to saturation 

input = 3 Sin 5.2irt, K = 70 
* * m 



83 




84 



FIGURE C-6: Fundamental oscillation to 1/3 subharmonic oscillation for relay. 

Input = 0.27 Sin 2irt 
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FIGURE C-7: 1/3 to 1/4 subharmonic oscillation due to relay with hysteresis, 

a = 1, c = 0.3; input = 0.2 Sin 1.87Tt to 0.2 Sin 2iTt. 




FIGURE C-8: 1/3 to 1/5 order subharmonic for system with 

relay; a = 1, c = 0.1. Input frequency from 
1.8 hz to 1.9 hz. 
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FIGURE C-9 : Change of input frequency from l.Shz to 1.3hz to 1.2hz resulting in a change 

of subharmonic oscillation from 1/5 to 1/4 to 1/3. a = 1, c = 0.1. 
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FIGURE C-10: 1/5 subharmonic oscillation to fundamental 

oscillation. a = 1, c =.01. Input fre- 
quency from 4.0 hz to 4.6hz. 
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FIGURE C-ll: 1/5 to 1/7 subharmonic oscillation for relay 

a = 1, c = 0.1. Input frequency from 3.1hz 
to 3.2hz. 
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FIGURE C-12: 1/7 to 1/5 subharmonic oscillation. Input 

frequency from 2.5hz to 2.4hz, a = 1, c = 0.1. 
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FIGURE C-13: Change of input frequency from 3.2hz to 3.4hz 

causing change of subharmonic oscillation from 
1/7 to 1/9; a = 1, c = 0.1. 
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FIGURE C-14: 1/9 to 1/7 subharmonic for relay, a = 1, c = 0.1. 

Input frequency from 3.5hz to 3.4 hz. 
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FIGURE C-15: 1/11 to 1/9 subharmonic oscillation. Input 

frequency from 4.3hz to 4.2hz. a = 1, c = 0.1. 
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FIGURE C-16: Output oscillating at two different subharmonics 1/11 and 1/13 

at same input frequency of 5.5hz. a = 1, c = 0.1. 




FIGURE C - 1 7 : Signal waveform of system with two nonlinearities 

C = 0.387. 

(a) Input = 9Sin 3.44iTt (d) Output of N, f(t) 

(b) Error Signal, E ( e ) i npu t to N ? , e. 

(c) Input to N l5 e (f) System Output 
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FIGURE C-18: Singal waveforms for 1/5 

lation. Input = 95in44 

(a) Input signal 

(b) Error signal 

(c) Input to first nonlinearity 



subharmonic oscil- 
t; K m =50, k t =0.04. 

(d) Input to second 
nonlinearity 

(e) System output 
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saturation and backlash. 
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FIGURE C-20: 1/3 subharmonic oscillation build up with input = 9Sin6.6-rrt for saturation 

and backlash, K, = 0.06, K = 200 
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FIGURE C-21: Build up of 1/7 subharmonic saturation and backlash. K. = 0.02, K = 100. 




FIGURE C-22: 1/3 subharmonic oscillation fading to funda- 

mental oscillation due to change of input 
frequency from 1.3 to 1.4 hz. Saturation 
and backlash: K = 0.13, K = 100. 
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FIGURE C-23: Signal waveform showing change of subharmonic 

mode from 1/4 to 1/5. 



(a) 
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= 21Sin3 . 2irt 
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